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Abstract 

Within the non-equilibrium Green's function technique on the real time contour, 
the ^-functional method of Baym is reviewed and generalized to arbitrary non- 
equilibrium many-particle systems. The scheme may be closed at any desired order 
in the number of loops or vertices of the generating functional. It defines effective 
theories, which provide a closed set of coupled classical field and Dyson equations, 
which are self-consistent, conserving and thermodynamically consistent. The ap- 
proach permits to include unstable particles and therefore unifies the description 
of resonances with all other particles, which obtain a mass width by collisions, de- 
cays or creation processes in dense matter. The inclusion of classical fields enables 
the treatment of soft modes and phase instabilities. The method can be taken as a 
starting point for adequate and consistent quantum improvements of the in-medium 
rates in transport theories. 



1 Introduction 



Non-equilibrium Green's function technique, developed by Schwinger, Kadanoff, Baym 
and Keldysh [1-4] , is the appropriate concept to study the space-time evolution of many- 
particle quantum systems. This formalism finds now applications in various fields, such as 
quantum chromodynamics [5], nuclear physics [6-12], astrophysics [10,13,14], cosmology 
[15], spin systems [16], lasers [17], physics of plasma [18,19], physics of liquid 3 He [20], 
critical phenomena, quenched random systems and disordered systems [21], normal metals 
and super-conductors [13,22,23], semiconductors [24], tunneling and secondary emission 
[25], etc. 
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For actual calculations certain approximation steps are necessary. In many cases pertur- 
bative approaches are insufficient, as for physical systems with strong couplings, e.g. like 
those treated in nuclear physics, in physics of liquids 3 He and 4 He, or high-temperature 
super-conductivity, etc. In such cases, one has to re-sum certain sub-series of diagrams 
in order to obtain a reasonable approximation scheme. In contrast to perturbation the- 
ory for such re-summations one frequently encounters the complication that the resulting 
equations of motion, even though self-consistent, may no longer comply with the con- 
servation laws, e.g., of charge, energy and momentum. This problem has first been con- 
sidered in two pioneering papers by Baym and Kadanoff [26,27] discussing the response 
to an external perturbation of quantum systems in thermodynamic equilibrium. Baym, 
in particular, showed [27] that any approximation, in order to be conserving, must be 
so-called ^-derivable. Thereby, he exploited the properties of an auxiliary functional, the 
^-functional, introduced by Luttinger and Ward [28] a year earlier for the formulation of 
the thermodynamic potential (see also [29]). Thereby the ^-functional is determined in 
terms of full, i.e. re-summed, Green's functions and free vertices. The scaling parameter of 
the vertices can be considered as an expansion parameter of a given approximation level. 
In the non-equilibrium formalism the problem of conserving approximations could be even 
more severe than in the case of the systems response to an external perturbation close to 
thermal equilibrium, since the system may exercise a rather violent evolution. Apart from 
transport models, mostly based on the quasi-particle approximation like Landau's Fermi 
liquid theory, there were only few attempts to discuss the issue of conserving approxima- 
tions in the context of the non-equilibrium field theory (see, e.g., [2,6,9]), which mainly 
considered Hatree-Fock and T-matrix approximations. However, the general problem of 
constructing conserving approximations in the non-equilibrium case and, in particular, 
beyond the quasi-particle limit has not explicitly been addressed yet. 

Alongside, the question of thermodynamic consistency is vital. If, as a result of a non- 
equilibrium evolution, a system arrives at an equilibrium state, the non-equilibrium Green's 
functions should properly describe thermodynamic quantities and potentials, such that 
thermodynamic relations between them are preserved. This problem is also relevant to the 
thermodynamic Green's function technique, as already considered by Baym [27]. Baym 
demonstrated that any ^-derivable approximation is at the same time thermodynamically 
consistent. 

In this paper we re-address the above problems and extend the concept i) to the genuine 
non-equilibrium case formulated on the closed real-time contour, ii) to relativistic field- 
theory Lagrangians in principle of arbitrary type (not just two-body fermion interactions) 
and iii) to the inclusion of classical fields, i.e. non-vanishing expectation values of the field 
operators. The generalized scheme permits to construct self-consistent, approximate, cou- 
pled dynamical equations of motion for the classical fields and Green's functions of the 
system on the closed real-time contour. This set of equations is conserving and thermody- 
namically consistent. Thereby, the inclusion of classical fields permits to account for the 
phase-transition phenomena or to describe the coherent dynamics of soft modes, much in 
the spirit of hard-thermal-loop re-summations [30,31,12]. Avoiding the quasi-particle limit 
enables us to appropriately consider the finite mass-width of all constituents in the dense 
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matter environment. The latter aspect unifies the description of resonances, which have 
already a mass-width in vacuum, with all particles, which acquire a dynamical width dur- 
ing the collision processes in the dense matter. The proper account of the finite width of 
the particles is also vital for a non-singular treatment of the soft-mode problem [12], where 
space-time coherence effects, like the Landau-Pomeranchuk-Migdal effect, defer the use 
of zero width quasi-particles and require non-perturbative re-summations. In this paper 
we confine the presentation to the derivation of the closed set of self-consistent couple 
Kadanoff-Baym and classical field equations. This constitutes the basis for various fur- 
ther steps towards classical-type transport schemes through the gradient approximation, 
which will be presented in a forthcoming paper [38]. 

For the sake of clarity, we restrict the presentation to systems of relativistic scalar bosons. 
This allows us to formulate the basic ideas in a simple and transparent form. The re- 
sulting relations can directly be generalized to multi-component systems of relativistic 
bosons and fermions. In sect. 2 we introduce the general equations of motion and the 
expressions for the conserved quantities on the operator level. The equations of motion of 
the corresponding expectation values are formulated within the real-time closed contour 
formalism (sect. 3). Thereby, it is advantage to formulate the concepts in terms of gen- 
erating functionals on the non-equilibrium contour, where the special functional $ plays 
a central role (sect. 4). This generating functional takes the same status in the space of 
Green's functions (two-point functions) and classical fields (one-point functions), as the 
original Lagrangian for the field operators. Subsequently, we formulate the diagrammatic 
representation for <3> (sect. 5). We show that any approximation, where all classical field 
sources and self-energies are ^-derivable in the sense of a variational principle, has the 
following properties: (i) it is conserving, (ii) it provides conserved current and energy- 
momentum tensor, which are identical to the corresponding Noether quantities (sect. 6), 
and (iii) it is at the same time thermodynamically consistent (sect. 7). In the summary, we 
formulate the main results and briefly discuss extensions and applications of the derived 
formalism. The list of diagrammatic rules is deferred to the Appendix A, while Appendix 
B contains some helpful equilibrium relations. 



2 Energy-Momentum Tensor and Conserved Currents 

We consider a system of relativistic scalar bosons, specified by the free Klein-Gordon 
Lagrangians 




for neutral bosons, 



for charged bosons, 



(2.1) 



where (f>(x) and (f>^(x) are bosonic field operators. The convention of units is such that 
h — c — 1. The interaction Lagrangians £ int {0} (for neutral bosons) and £ mt {0, <fi^} 
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(for charged bosons) are assumed to be local, i.e. without derivative coupling. Under these 
conditions the Lagrangians are charge symmetric. Charges are understood as the electric 
charge, strangeness, iso-spin, etc. 

The variational principle of stationary action leads to the Euler-Lagrange equations of 
motion for the field operators 



-S x <j>(x)=J(x) = -^, where S x = -d^-m 2 , (2.2) 

and similarly for the corresponding adjoint equation. Thereby, the J(x) operator is a local 
source current of the field 6, while S x is the differential operator of the free evolution with 
the free propagator A°(y,x) as resolvent. 

The standard canonical energy-momentum tensor [32] has some undesired features, as 
it is non-symmetric in the Lorentz indices, for example. Alternatively, using the Euler- 
Lagrange equations of motion and the definition of the source current (2.2), one can show 
that the following form also defines a conserving energy momentum tensor 



U (alt.) 1 



V 

2 / neut. 



(5 - d[ 



+g »v (2 in \x) - £ pot (x 



with d^ lt) (x) 



0. 



dC°(x) 



6{y) - 6\x) 



x=y 



(2.3) 



For notational simplicity, expression (2.3) and similar expressions below, which appear 
symmetric in 6 and 6\ are written in such a way that they directly apply to complex 
fields. The symbol (l/2) neut implies that for real fields the corresponding expressions are 
obtained from those for complex fields by multiplying by the factor 1/2 upon equating 
= 6. Above, we have introduced operators of the interaction energy density £ mt (x) of 
the system, which accounts for the total interaction part of the energy density, and the 
potential energy density £ pot (x), both given as 



£ mt (x) = -£ mt (x), 

^'nGLl^^^) (2 - 4) 

=-\(\) nM {jHx)i(x)+mx)j( X )). 

For a multi-component system, the latter defines the sum of the potential energy densities 
of any field 4>{x) with the currents J(x) induced by the other fields in the system. 
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In terms of the differential operator p% = id%, the alternative energy-momentum tensor 
of eq. (2.3) can be written in a charge symmetric form as 



1 /l 



e (Tit.)(^-4 ^ 



neut. 



x=y 



+ g^ u (S int (x) -S pot (x)) . (2.5) 

This form of the energy-momentum tensor is equal to the metric one, which results from 
variation of the action over the metric tensor g^ u rather than over the fields. All terms in 
eq. (2.5) are evidently symmetric in \iv. 

For specific interactions, eq. (2.4) provides simple relations between S mt (x) and £ pot (x). 
If all vertices of £ mt have the same number 7 of field operators attached, one simply 
deduces 



t n \x) = -£ pot (x). (2.6) 

7 

For instance, for the 4 -theory, where 7 = 4, the interaction energy is half of the potential 
energy. 

If the Lagrangian is invariant under some global transformation of charged fields (with 
the charge e), e.g., 

$( x ) e~ ieA 0(x); $\x) e ieA 0t( x ) ) ( 2 .7) 
there exists a Noether current defined as [32] 

(2.8) 



dC(A(x),d,A(x)) 

J (Noether)- 8 (8^A{x)) 



A(x)=0 



which is conserved, i.e. <9 M j (* Noether ) = 0. Formally, it is derived by applying the local 
transformation of the form (2.7) and using the stationary condition of the action around 
physical solutions. 

From the Euler-Lagrange equations of motion (2.2), one obtains 

^Noether) = (^tofts) - ^ (x) J (x)) , (2.9) 



which vanishes, since for the symmetry (2.7) the interaction Lagrangian has to consist 
of terms, where <f>(x) and 4>\ x ) appear pairwise. In terms of the = \d% operator, the 



5 



i(Noether) curren t can again be written in a charge symmetric form as 

^Noether)^) = ^ [((%)*+ &£) (P(x)$(y) + $ (y) & (x))] ^ , (2.10) 



which naturally vanishes for the neutral particles (e = 0). 

One may also define the tensor M^ vp , which is associated with the Lorentz invariance of 
the Lagrangian and provides the angular momentum conservation. However, we do not 
treat this tensor in this paper, since it is of no common use in kinetics. 



3 Real-Time Contours 



In the non-equilibrium case, one assumes that the system has been prepared at some 
initial time to described in terms of a given density operator p — Ha Pa \o>) (o>\, where 
the | a) form a complete set of eigenstates of pQ. All observables can be expressed through 
n-point Wightman functions of Heisenberg operators A (ti), . . . , O (t n ) at some later times 



(6(t n ) . . . BfaAfa)) =: Tr 0(t n ) . . . B(t 2 ) A(ti)p (*o) 

= ]T P a (a\ 0(t n ) . . .B(t 2 )A(t!) \a) . (3.1) 

a 

Note the fixed operator ordering for Wightman functions. 



t tt 

Figure 1: Closed real-time contour with two external points x,y on the contour. 

The non-equilibrium theory can entirely be formulated on a special contour — the closed 
real-time contour (see figure 1) with the time argument running from t to oo along the 
time- ordered branch and back to to along the anti-time- ordered branch. Contour-ordered 
multi-point functions are defined as expectation values of contour ordered products of 
operators 

(T c A(x 1 )B(x 2 ) ...) = (tcM^B^) . . .exp |i| 2j nt dx , (3.2) 
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where Tc denotes the special time-ordering operator, which orders the operators according 
to a time parameter running along the time contour C. The l.h.s. of eq. (3.2) is written in 
the Heisenberg representation, whereas the r.h.s. is given in the interaction (I) representa- 
tion. Here and below, the subscript " I" indicates the interaction picture. Functions with n 
points can be expressed in terms of products of other multi-point functions contour inte- 
grated over internal coordinates. Ultimately, one likes to express the multi-point functions 
of interest in terms of one- and two-point functions (Wick's linked cluster expansion), i.e. 
in terms of classical fields, classical source-currents, propagators and self-energies. Note 
that at this level the contour is not a contour in the complex plane, as the figure may sug- 
gest, but rather it runs along real time arguments. It is through the placement of external 
points on the contour that the contour ordering obtains its particular sense. 

In certain calculations, e.g., in those that apply the Fourier and Wigner transformations, 
it is necessary to decompose the full contour into its two branches — the time- ordered and 
anti-time- ordered branches. One then has to distinguish between the physical space-time 
coordinates x, . . . and the corresponding contour coordinates x c which for a given x take 
two values x~ = (x~) and x + = (/i G {0, 1, 2, 3}) on the time ordered and anti-time 
ordered branches, respectively (see figure 1). Closed real-time contour integrations can 
then be decomposed as 

oo to oo oo 

Jdx c ... = Jdx~ ...+ J dx + ... = J dx~ . . . - J dx + . . . , (3.3) 

C t oo t t 

where only the time limits are explicitly given. Thus, the anti-time-ordered branch acquires 
an extra minus sign if integrated over physical times. For any two-point function F, the 
contour values are defined as 

F i i(x,y):=F(x i ,y j ), i,j e {-,+}, (3.4) 



on the different branches of the contour. The contour 5-function is defined as 



8 c (x % ,y 3 ) = <T %3 5*(x-y), a 1 ' 



1 
-1 



(3.5) 



where, <7ifc is the 03-Pauli matrix. For any multi-point function, the external point x max , 
which has the largest physical time, can be placed on either branch of the contour without 
changing the value, since the contour-time evolution from x^ nax to provides unity. 
Therefore, one-point functions have the same value on both sides of the contour. 

Due to the change of operator ordering, genuine multi-point functions are discontinuous in 
general, when two contour coordinates become identical. In particular, two-point functions 
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like iF(x,y) = (T c A(x)B(y)^ become^] 



iF(x,y) 



iF (x,y) iF + (x,y)^ 
\iF+-(x,y) iF ++ (x,y) J 



/ (TA(x)B(y)) (B(y)A(x)) > 

v (A(x)B(y)) (T^A(x)B(y))j 



(3.6) 



where T and T _1 are the usual time and anti-time ordering operators. Since there are 
altogether only two possible orderings of the two operators, in fact given by the Wightman 
functions F h and F + ~, which are both continuous, not all four components of F are 
independent. From eq. (3.6) follow relations between non-equilibrium and usual retarded 
and advanced functions 



F R {x,y) = F—(x,y) - F- + (x,y) = F + ~(x,y) - F ++ (x,y) 
:= 0(^o - yo) (F + ~{x, y) - F~+{x, y)) , 

F A (x,y) = F-(x,y) - F+-(x,y) = F- + (x,y) - F ++ (x,y) 
:= -e{y - x ) (F + -(x, y) - F-+(x, y)) , 

where 0(xq — yo) is the step function of the time difference. 



(3.7) 



Discontinuities of a two-point function may cause problems for differentiations, in particu- 
lar, since they often occur simultaneously in products of two or more two-point functions. 
The proper procedure is, first, with the help of eq. (3.7) to represent the discontinuous 
parts in F and F ++ by the continuous F~ + and F + ~ times G-functions, then to com- 
bine all discontinuities, e.g. with respect to xo — yo, into a single term proportional to 
Q(xo — yo), and finally to apply the differentiations. One can easily check that in the 
following particularly relevant cases 



/ dz (F{x\z)G{z,x i )-G{x i , z)F(z, x J 



_d_ 

dx, 



J dz (F(x\ z)G(z, x j ) + G(x\ z)F(z } x j 



d d 



dxa dy u 



dz (F(x\z)G(z,y j )-G(x i ,z)F(z,y j ) 



(3.8) 
(3.9) 

(3.10) 



x=y 



all discontinuities exactly cancel. Thereby, these values are independent of the placement 

[ F c F<\ 

Quite commonly, like in refs. [2,6], the notation F = I I is used for two-point func- 

\p> pa J 

tions instead of (3.6). We prefer the more flexible { — h} labeling of contour points. 



8 



of x l and on the contour, i.e. the values are only a function of the physical coordinate 



x. 



Boson fields may take non- vanishing expectation values of the field operators (f)(x) = (</>V 
called classical fields. The corresponding equations of motion are provided by the ensemble 
averaging of the operator equations of motion (2.2) 

S x <j>(x) = -J(x), or <f>(x) = <t>°{x) - J dyA°(x,y)J(y). (3.11) 



Here J(x) = (^ J(x)j, while <p°(x) = i^(f)i{x)j is the freely evolving classical field which 
starts from (j)(to,x) at time t - Thereby, A°(x,y) is the free contour Green's function 

\A°(x,y) = (Tcfaxtftiy)) - <f>°(x)(<f>°(y)y. (3.12) 



The reader can easily verify that the contour form (3.11) is equivalent to the standard 
retarded form of the classical field equations due to eq. (3.7) and the fact that J(x) and 
<f>(x) are one-point functions, which have identical values on both sides of the contour. 
The free propagator is resolvent of equation 

S x A°(x,y) = S c (x,y) (3.13) 

on the contour, where Sc(x\ y- 7 ) is the contour 5-function (3.5). Subtracting the classical 
fields via 

$ = <P+$, (3.14) 
we define the full propagator in terms of quantum-fluctuating parts (p of the fields 



iA(x,y) = (T c ip{x)ip\y)) = (T c $(x)$*(y)) - <t>{x)<t>*{y) = (T C $(x)$\yj) e . (3.15) 

Here and below, the sub-label "c" indicates that uncorrelated parts are subtracted. In 
terms of diagrams it implies, that the corresponding expectation values are given by sums 
of entirely connected diagrams. 

Averaging the operator equations of motion (2.2) multiplied by ffl(x) and subtracting 
classical-field parts, one obtains the equation of motion for the propagators as 

S x A(x, y) = 5 c (x, y) + i (T c J(x) <p\y)) c , (3.16) 

where the contour 5-function appears due to the contour ordering of the operators in 
A(x,y). 
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Eq. (3.16) is still exact and accounts for the full set of initial correlations contained in 
Pq. In order to proceed, one may suggest that the typical interaction time r int for the 
change of the correlation functions is significantly shorter than the typical relaxation time 
T re i, which determines the system evolution, when one neglects those initial correlations. 
Then, describing the system at times t — to 3> Tint, one can neglect the initial correlations 
which are supposed to be dying at the time ~ r int in accordance with the Bogolyubov's 
principle of the weakening of initial correlations.^ As a result one can apply the standard 
Wick decomposition dropping higher order correlations for the driving terms on the r.h.s. 
of both equations of motion (3.11) and (3.16). Then both driving terms can be expressed 
solely as functionals of the classical fields and one-particle propagators rather than on 
higher order correlations. In particular the driving term of (3.16) can then be expressed 
through the proper self-energy II 



i(T c J{x)(p\y)) =i(T c 




{'/ 


dz£j nt 1 








d 




exp | 


d<pi{ 


z) 



J dzU(x,z)A(z,y). 



J f> int 



T c ${z){p\y) 



cl 



(3.17) 



The second line of eq. (3.17) results from the fact that the expectation value is connected 
and, therefore, in the first line, <p\(y) has to be contracted with one of the operators 
(pi(z) occurring in exp ji J c dz'£J nt | J\(x). The expression for the self-energy II deserves 
attention, since we have separated the full propagator in (3.17). In order to achieve the 
proper counting of terms, II has to be one-particle irreducible (label cl), i.e. the corre- 
sponding diagram cannot be split into two pieces, which separate x from z by cutting 
a single propagator line. Obviously, II may have singular (5-functional) one-point parts 
and genuine two-point parts (the latter are given by all connected one-particle irreducible 
diagrams of the current-current correlator) 



-m(x,y) = {T c / Sc(x,y)-(T c J(x)j\y)) ic , (3.18) 

in the Heisenberg picture. 

With the help of eqs (3.16) and (3.17), one recovers the Dyson equation in the differential 
form 



2 Actually, considering a dilute gas limit, Bogolyubov suggested the weakening of all the cor- 
relations, whereas we use a weaker assumption on the weakening of only short-range (~ T[ nt ) 
correlations, cf. [33]. 
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S x A(x,y) = 5c(x,y)+ J dzU(x, z)A(z,y), (3.19) 

c 

(S y )* A(x, y) = 5 c (x, y) + J dzA(x, z)U(z, y). (3.20) 

c 

Using the resolvent property of the free propagator, we can write down these equations 
also in the integral form. 

In diagrams, free and full classical fields are represented by "pins" with cross and "o-cross" 
as heads 



<f>°(x) = f 4>{x) = f , (3.21) 
while free and full propagators are given by thin and thick lines, respectively, 



y) = iA(x, y) = — - (3.22) 

Thereby, complex fields carry a sense, the arrow always pointing towards the <f> in the 
contour ordered expressions. In diagrammatic representation, the classical-field equations 
(3.11) and Dyson equations (3.19) are given by 




iJ 

f + Y , (3.23) 




+ ^ — (-ill) — «— (3.24) 



with the one- and two-point functions iJ(x) and — iU(x, y) as driving terms. 

The averaged values of conserved quantities can be expressed in terms of the one- and two- 
point functions introduced so far. Averaging the operator value of the energy-momentum 
tensor of eq. (2.5), we arrive at 



ncut. 



{m*+Py) (m*+%) U*(x)<f>(y)+iA s ^(y,x 



x=y 



+ (S' m \x) - S pot (x)) , (3.25) 
where S mt (x) = — (C mt (x)\ and the potential energy density becomes 
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S^(x) = (£v°\x)) = \ (i)^ | - \r{x)<P{x) + J(x)<r(x)] 

+ ijdz [U(x, z)A(z, x) + A(x, z)U(z, x)] j . (3.26) 

due to eqs (2.4) and (3.17). Note that we do not prescribe any contour indices to x in the 
integral term of £ pot , since actually this term is independent of the contour placement of 
x due to discontinuity property (3.9). The Noether current (2.10) takes the form 



MX) =■ (^Noether) H +3?) U* (^(y) + iA^(y , x) 



(3.27) 



x=y 



In order to keep the expressions and charge symmetric we have introduced the 
symmetric quantities 

F^(x, y) = \ (F(x~,y + ) + F(x + , „")) , (3.28) 



where F(x, y) is any two-point function on the real-time contour. This is not automat- 
ically provided by the variational methods leading to the Noether energy-momentum 
tensor and current, as they only provide the vanishing of the corresponding divergence. 
The integrated form of the conserved quantities has to be adjusted such that charge sym- 
metry is maintained, thus describing contributions of both, particles and anti-particles, on 
equal footing. This way one properly accounts for the modification of the vacuum polar- 
ization in the medium, since the vacuum-polarization energy coincides with the zero point 
energies of the field oscillations. The corresponding divergence has still to be appropriately 
renormalized. 



4 Functional W, T and $ 



The standard generating functional for connected n-point functions is given by the log- 
arithm of the expectation value of the time-evolution operator with external one-point 
sources r)(x) and rj* (x) 



W{r],r]*} = -iln(Texp 



= — i In ( exp 



ifdx 
T exp 



r)(x)4>\x) + rf(x)4>{x) 



i J dtH\ Texp i J dx(Z\ lt + [rj(x)$\(x)+r]*(x)$ 1 (x)\) 



where H j is the free Hamiltonian of the system in the interaction representation. Since 
the functional dependence concerns only the external sources rj(x), the operator part can 
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be cast into different pictures, such as the Heisenberg or interaction ones. The latter 
establishes the perturbation expansion. Here and below, (...) denotes a trace over all 
states, which includes the ensemble average over the density operator p — p(t ) at 
initial time t (cf. eq. (3.1)). 

It is advantageous to introduce a scale function X(x), which scales interaction vertices 
of the interaction Lagrangian density £ int at space-time coordinate x and defines a X(x) 
dependent interaction Lagrangian density 

CT = X(x)C^{^(x)J(x)}. (4.1) 



This scaling provides the clue to the proof of the diagrammatic representation of the 
auxiliary functional $ in terms of closed diagrams. 

We generalize the above W functional to the real-time contour C. Following Luttinger and 
Ward [28] and Baym [27], we also extend it to include external bilinear sources K(x,y), 
besides the interaction scale X(x), 



W {i], K, X} = -i In ^exp 



-i / dxH^ 
c 



T c exp<(i J dx(XC\ n 

c 



V 

v ^ / neut. 



r](x)$\(x) +r}*(x)$i(x) +i J dy$i(x)K(x : y)$\(y) 



(4.2) 



While n-th order functional variations of (4.2) with respect to 77 generate the connected 
n-point functions, the first functional variation with respect to K(x, y) gives the total 
two-point propagator including disconnected pieces 



5W{ V ,K,X}=[- 

v ^ / neut. 



dx [(f)*(x)5r](x) + (f)(x)Srj*(x)] 



+ jdxjdy [<j>(x)<j>*(y) + A(x,y)]6K(y,x) 



c c 

int / 



-J dx£™(x)5X(x)/X(x). 
c 



(4.3) 



Here, and A denote the classical fields and propagators, respectively, which now implic- 
itly depend on A. The term, resulting from the variation of the interaction scale function 
X(x), defines a one-point function S 1 ^ t (x) = — ^^(x)^, which agrees with the corre- 
sponding expectation value of (2.4) but for the scaled Lagrangian. 



The step towards a functional that depends on the classical fields and propagators 
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A rather than on the external sources 77 and K is provided by the double Legendre 
transformation of W to r{0, 0*, A, A} given as 



T{ ( f ) , ( f ) *,A,X} = W{r ] ,K,X} 



dx [rj*(x)(j)(x) + r/(x)(j)*(x)} 



+ / dxdy [(f)(y)(f)*(x) + A(y,x)]K(y,x) 
c 



(4.4) 



Here, the sources 77 and K have to be expressed through and A. Apart from the 5 A 
dependences, the functional variation 



5r{0,0*,A,A} = -(i 

k ^ / neut. 



dx \rf(x)8<f)(x) + rj(x)S<p*(x)] 



+ J dxj dy [8<f>{y)4>*{x) + (j)(y)8(j)*(x) + 6A(y, x)] K(y, x) 
c c 

dxS[ nt {x)5X(x)/X{x) 



(4.5) 



vanishes at vanishing external sources 77 and K. The latter together with the condition 
A = 1 corresponds to the physical solution. Note also that for the physical solution, i.e. at 
i] = K = 0, the values of the two functionals T and W are identical. Indeed, the variations 



ST/54> = 0, 5T/54>* = 0, 5T/5A = 



(4.6) 



provide us with equations of motion for classical fields (3.11) and the complex conjugated 
one, as well as Dyson's equation (3.19). 

Following Luttinger, Ward [28], and later by Cornwell et al. [34], who used path-integral 
methods for case of the imaginary-time formulation of equilibrium systems, we represent 
our T functional related to the real-time quantities in the form 



T{0, 0*, A, A} = T° + J dx£ {<P, «9 M 0} 



In (l - ©A n) + ©A n] + $ {0, 0*, A, A} 



+1( 9 

v ^ / neut. 



(4.7) 



this way defining the auxiliary functional $ {0, 0*, A, A}[^} Constructing $ solely in terms 
of the fields and one-particle propagators complies with the assumption of ignoring higher 



3 Please note that our $ is different from the auxiliary quantity T2 defined by Cornwell et. 
al. [34] to the extent that their is void of zero- and one-loop terms. These terms are rather 
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order correlations. The T and C° parts, where C° is the free classical Lagrangian function, 
represent the non-interacting parts of T. Thereby, the T term solely depends on the 
unperturbed propagator A and hence is treated as a constant with respect to functional 
variations of T. The ln(. . .) is understood in the functional sense, i.e. by a series of re- 
folded contour convolutions, denoted by the ©-symbol, formally resulting from the Taylor 
expansion of the ln(l + x) at x — 0. This ln-term accounts for the change of T due to the 
self-energies of the particles. The T , C° and In terms in eq. (4.7) account for the one-body 
components in the T. The remaining A II and $ terms correct for the true interaction 
energy part of T. As shown in the next section, the form (4.7) presents a re-summation 
of the corresponding perturbative expansion of the value of T (cf. eq. (4.4)) in terms of 
full classical fields and propagators. 

The specific form (4.7) has important functional properties, which provide us with a 
number of useful relations. Functional variation of T {(f), <f>*, A, A} in the form of eq. (4.7) 
leads to 



5T {(f), 0*, A, A} = ( - ) I I dx [5(f)(x)(S x )*(j)*(x) + 5(f)*(x)S x (f)(x) 

- / neut. | J 

- l — — - © A° - ©a) © 511 

1 - ©A IT J 



_!(,. 

+ i 



J dxdyU{x, y)5A(y, x) \ + 5$ {0, 0*, A, A} . (4.8) 



Here 511 is understood as a variation induced by 6 A, 5(f), 5(f)* , and 5X, respectively. Eqs 
(4.5) and (4.6) imply the following variational rules for the auxiliary $ functional 



5 $ 



{>,0*,A,A}= Q)^ t Udx[J*(x)5(f)(x) + J(x)5(f)*(x)} 



- i j dxdyU(x,y)5A(y,x) \ - J dxS mt (x)5\(x), 
c ) c 



(4.9) 



or 



U(x) 
iU(x,y) 



5i$ 

5(f)* (x) ' 

5i$ 

5iA(y,x) 



x < 



2 for neutral bosons, 
1 for charged bosons, 



(4.10) 
(4.11) 



included in defining a full classical Lagrangian that depends on the classical fields and tad-poles, 
while we have placed all interaction parts into <5. 
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(4.12) 



The virtue of the functional form (4.7) is that these requirements can be met simultane- 
ously and that there exists a unique form of $, for which the three derived quantities — the 
one-body source current J(x), the two-body self-energy —iH(x,y) and the interaction en- 
ergy density S mt (x) — take their physical values at the physical solutions of the equations 
of motion (3.19) and (3.11). This will become clear in more detail in the next section, 
where we discuss the diagrams defining the various functionals. In its turn, $ can be seen 
as a generating functional for the source terms J of classical fields and self-energies II for 
the set of Dyson equations. Therefore, approximation schemes can be defined through a 
particular approximation to $. Thereby, the invariance properties of $ play a central role 
to define conservation laws for the approximate dynamics. 

It is important to emphasize that we do all functional variations independently of any place 
on the contour. Thus, different contour times are considered as independent, even though 
they may refer to the same physical timej^]. In principle, all variational considerations 
given in this section apply to any kind of time contour, even to non closed and complex 
ones as well as to any operator p defining the averaging (...) = Tr {.. . p } including the 
unit operator, as used in Matsubara's imaginary-time formalism. For a particular choice 
of p and of a contour the physical values of W and T are identical for the corresponding 
physical solutions along this contour. In the imaginary-time method the value V = W 
takes the meaning of the thermodynamic partition sum. In this paper we concentrate on 
the non-equilibrium closed real-time formalism for which Trp = 1, and therefore the 
physical values of W and Y trivially vanish, i.e. 

W = T = $ = for physical solutions of A, 0, 0* on the contour. 

An important comment must be given at this stage. One should clearly distinguish be- 
tween the functional form of a functional, which acquires its meaning through variational 
methods, and the physical value that functional takes once the physical solutions of the 
equations of motion are inserted. For instance, two functionals W and T are completely 
different in their functional meaning, while they take the same physical value for the 
physical solution. Therefore, for all functionals the functional dependences are explicitly 
given in braces. Our strategy below will be first to perform general variations of T and 
$, allowing non-physical values of A, <p, <ft* and A, and only then to put them to their 
physical values. This way, a number of important relations between Green's functions, 
self-energies and mean fields will be obtained. 



The fact that for the physical solutions the components of A on the different branches of 
the contour are not independent (cf. (3.7)), has no importance for the variational procedure. 
The reason is that rules (3.7) only apply to the physical A and <fi, which are provided by the 
stationary "points" of the variational principle, i.e. solving the equations of motion (3.11), (3.19) 
and (3.20). 
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5 Diagrams for T, $ and S™ l {x) 



According to eqs (4.5) and (4.12), we have 



dxS 



int / 



[X 



A— r{0{A},<f{A},A{A},A} 



A=l 



A^-${0,0*,A,A} 



(5-1) 



J A=l 



where now A is treated as a global scale parameter (note that only a partial derivative is 
applied to $, i.e. the <p, (f>* and A values are kept constants.). In the perturbation theory, 
the diagrammatic rules to calculate the one-point function S mt (x) are straightforward 




Here the diagram symbolically denotes all connected (label c) closed perturbation-theory 
diagrams generated by expanding the exponential function in (5.2). The full dot denotes 
the external point x which is not integrated out. Integrating (5.1) with respect to A, we 
define the quantity iT given by the following perturbative diagrammatic representation 




where the integration constants have been chosen such that for physical solutions F = T. 
One can see that each diagram contributing to f has to be weighted with its inverse 
number of vertices l/n\, due to the formal A-integration of (5.2). It is important to realize 
that due to these global factors such a set of diagrams is not resumable in the standard 
diagrammatic sense^]. Also T in the form of eq. (5.3) is a functional of 0°, 0°*, A , A rather 
than of (f>, (ft*, A, A, as needed for T and $ as discussed in the previous section. However, 
we can arrive at the required functional dependence of T as follows. The expression (5.2) 
for — iS mt (x) can be re-summed and entirely expressed in terms of full classical fields and 
full propagators. The re-summed diagrams are then void of any self-energy insertions and 
therefore have to be two-particle irreducible 



-i£ mt (x)=£ 









) 












^c2 



5 Diagrammatic re-summation implies that sub-diagrams with the same external structure (i.e. 
same number of external points and types of propagators to be attached at each external point) 
can be summed up to give a total re-summed expression that can then be embedded into more 
complicated diagrams, e.g. self-energy insertions can be re-summed to full Green's functions. 



17 



Diagrams of class c2 cannot be decomposed into two pieces by cutting two propagator 
lines. The formal integration of the last equality in (5.1) with respect to A keeping and 
A constant provides the diagrammatic expression for $ in terms of full Green's functions 
and classical fields. Therefore, ir {0, 0*, A, A} can be expressed in terms of the following 
diagrams (cf. eq. (4.7)) 

ir {0, 0*, a, A} = ir° {a } + i J dxc ^, ^0} 

c 




- In (l - ©a° ©n) -oAon +i$ 



Here nn counts the number of II insertions in the ring diagrams providing the ln-terms, 
while for the closed diagrams of $ the value n x counts the number of vertices building 
up the functional $. Contrary to the perturbative diagrams of ir, cf. eq. (5.3), here 
the diagrams contributing to $ are given in terms of full propagators A and full time- 
dependent classical fields 0. As a consequence, these diagrams have to be two-particle 
irreducible (label c2). The latter property is required because of the re-summations of 
S mt (x). This also matches the diagrammatic rules for the re-summed self-energy U(x,y), 
which results from functional variation of $ with respect to any propagator A(y,x). In 
graphical terms, this variation is realized by opening a propagator line in all diagrams 
of $. The resulting set of thus opened diagrams must then be that of proper skeleton 
diagrams of II in terms of full propagators, i.e. void of any self-energy insertion. 

The diagrammatic rules for $, £ mt (x), J and II are determined by the following steps: 

(a) For all bosonic fields in i£ mt , replace by + (p in order to account for the classical 
fields (cf. (3.14)); 

(b) consider all possible pair contractions of the field operator (p(x) with <p^(y) in the 
formal expressions (5.6)-(5.9) given below and replace them by iA(x, y); 

(c) keep only those terms that correspond to two-particle irreducible diagrams for $ , i.e. 
which cannot be split into two pieces by cutting two different propagator lines. 

Further details are given in Appendix A. The diagrams of i$, — i£ mt (x) , iJ(x) and 
— in(x, y) are then generated by applying the above general rules to the following for- 
mal expressions 
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dx n (T c iC mt (x 1 )...iC mt (x n ) 



c 

> int / 



2c{A} ' 



-i£ int (:r) = ( T c iC mt {x) exp ( i J dx'C int {x') 



2c{A} 



(5.6) 



(5.7) 



u{x) = \ Tc W(x) ew VJ dxC {x) 

V C / ' 2c{A} 



(5.8) 



2c{A} 



where the sub- label 2c {A} refers to the above point (c). 



(5.9) 



As an example, we quote the diagrams in neutral scalar g0 4 /4! theory. Up to two vertices, 
the functional $ is given by the following expressions 



-15 



4! 



dx ((j) 4 (x) + Q(f) 2 (x) {(p{x)(p{x)) c + 3 (ip (x) (p (x))f) 



jdxjdy(A- 4!0(x)0(y) (^(x)^(y)>; + 4! (0(s)0(i/)>;) 



2 V 4! 

c c 

or in terms of diagrams (cf. also Appendix A) 



(5.10) 



+ ..., 



i$ 




+ 




' 1 " 




1 


.4!. 




_2-2! 



2 2 -2! 



" 1 " 




' 1 " 


.3!. 







3 " 
(5.11) 



The l/n\ factors are explicitly given, while the combinatorial factors according to rule (vii) 
in Appendix A are given in square brackets below each diagram. Functional derivatives 
with respect to <fi (pins) and propagators (full lines), cf. eq.(4.9), determine the source 
J(x) of the classical field and the self-energy U(x, y), respectively, 
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' 1 " 




r 




' 1 ' 


.3!. 




_2_ 




.3!. 



(5.12) 



-ill(x,y) = 




' 1 ' 




l" 




" 1 " 




' 1 " 


.2!. 




_2_ 




.2!. 




.3!. 



Small full dots define vertices which are to be integrated over, while big full dots specify 
the external points x or y; the first two diagrams of H(x, y) give the singular 8c (x, y) parts 
arising from classical fields and tad-poles. 



6 ^-Derivable Approximations and Invariances of $ 



The expressions for W, T and $ given so far are exact and represent a convenient formula- 
tion of the theory in terms of full propagators and self-energies. However, for any practical 
calculation one needs certain truncated approximate schemes. In the weak-coupling limit, 
one can restrict the perturbation series for the Green's function to a certain order. Then, 
as far as conservation laws are concerned, one encounters no particular problems, as con- 
served quantities are conserved order by order in perturbation theory. On the other hand, 
such perturbative expansion may not be adequate, as, for example, in the strong coupling 
limit, where re-summation concepts have to be applied, which re-sum certain sub-series of 
diagrams to any order. For such re-summation schemes, the situation with conservation 
laws is not that obvious. 

We consider so-called ^-derivable approximations, first introduced by Baym [27] within 
the imaginary time method. Such approximations are constructed by confining the infi- 
nite diagrammatic series for $ either to a set of a few diagrams or to some sub-series 
of diagrams. Note that the approximate $( a PP r ) itself is constructed in terms of "full" 
Green's functions and "full" classical fields, where "full" now implies that we have to 
self-consistently solve the classical-field and Dyson equations with the driving terms de- 
rived from this $( a PP r -) through relations (4.10) and (4.11). It means that even restricting 
ourselves to a single diagram in $( a PP r ), in fact, we deal with a whole sub-series of dia- 
grams in perturbation theory. Thereby, the term "full" takes the sense of the sum of this 
whole sub-series. Thus, a ^-derivable approximation offers a natural way of introducing 
closed, and therefore self-consistent approximation schemes based on summation of dia- 
grammatic sub-series. In order to preserve the symmetry of the exact $ with respect to 
permutations among i£ mt (xi) . . . i£ int (a; n ) (see eq. (5.6)), we postulate that a ^-derivable 



20 



approximation either takes into account all the diagrams up to a certain order n (in terms 
of "full" Green's functions and classical fields) or confines the treatment to only those 
diagrams with topologically equivalent vertices. As a consequence, approximate forms of 
$(appr.) d e fi ne effective theories, where $( appr ) serves as a generating functional for the 
approximate source currents j( appr -\x) and self-energies T[( appr -\x,y) (see eqs (4.10) and 
(4.H)) 



iJ (appr - } (x) 



£i$(appr.; 



S (0( a PP r -)*(x)) ' 



(6.1) 



^j<|)(appr.; 



5iA( appT -\y,x) 



x < 



2 for neutral fields, 
1 for charged fields, 



(6.2) 



which then are the driving terms for the equations of motion for the classical fields and 
propagators. The approximate $ also provides the corresponding expression for S mt (see 
eq. (4.12)). Below, we omit the superscript "appr.". 

We now like to demonstrate that ^-derivable approximations possess a number of remark- 
able properties. For such approximations, the invariances of <£> play as central a role as 
the invariances of the Lagrangian for the full theory. Thereby, the variational principle, 
where the interaction strength X(x), the classical fields <p(x), and propagators A(x,y) can 
be varied independently, provides a set of useful identities and relations. 

A general invariance of $ is provided by the substitution x x + £(x) for all integration 
variables in the contour integrations defining $. The Jakobi determinant required for each 
integration variable can be accommodated by a simultaneous change of the scale function 
X(x) at each vertex. Thus, the simultaneous variation 



<f)(x)^(f)(x + t(x)), 

A(x, y) A(x + £(x),y + £(y)), (6.3) 
A(x) = 1 =► X(x) = det ^ + , i.e. 5X(x) = 

leaves $ invariant. This way, one deduces 
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Interchanging x and y in the second II term in squared brackets, using partial integration 
and that the transformation £(x) can be chosen arbitrarily, one obtains the following 
relation 



dx 



— i 



2/ neut. 
^ / neut. 



x 



dx, 



dx^ j 



dx. 



dx,, 



(6.5) 



This is the key relation to prove energy-momentum conservation. It has features similar 
to a Ward identity, as it links derivatives of one-point functions with those of two-point 
functions. The two-point function contribution to this expression is of type of eq. (3.10), 
so that in eq. (6.5) the differentiations of the discontinuities, indeed, cancel out. 

With the help of the equations of motion (3.11), (3.19) and (3.20), the divergence of the 
kinetic term of the energy-momentum tensor (3.25) can be cast into 



1 n 



2 V2 



d„ 



neut. 



x=y 



-d tt <rE^(x)-(l) { [J(.rRo (,•) f JlMM.r) 

\£>J neut. 

+i J dz [U(x, z) ■ d*A{z, x) + &^A(x, z) ■ U(z, x)} 



(6.6) 



The remaining non-full-derivative terms on the r.h.s. of this equation are undesired, since 
they prevent us from presenting of the entire equation in the form of a full divergence. 
However, they can also be transformed into a full derivative form by means of identity 
(6.5). Adding this identity to eq. (6.6) cancels out the undesired last braced terms. Gath- 
ering all surviving terms of (6.6), we recognize it as the energy-momentum conservation 
law d^Q^lx) = with the energy-momentum tensor given by the Noether expression 
(3.25). Hence, the existence of a conserved energy-momentum tensor is proven for any 
^-derivable approximation. 

Along similar lines charge conservation can be proven, assuming that $ is invariant under 
the following simultaneous variation of classical fields and propagators 

<j,{x) e-' ieA{x) <f){x), <P*{x) e ieA{x) <f)*{x), A(x, y) e^ ieAW A(i, y)e ieA{y \ (6.7) 
Applying the rule (4.9) of the $ variation, to linear order in the phase A, one obtains 



5$ = eJ dx [J*(x)(f)(x) - J(x)(f)*(x)] [-iA(ar)] 
c 
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+ ie J dxdy [U(x, y)A{y, x) - A(x, y)U(y, x)} [-iA(x)] . (6.8) 
c 

Here we have used the variation rules for complex fields. Equating <5$ = for arbitrary 
A, we arrive at the identity 

e [J*(x)<j>(x) - J(x)<j>*{x)\ +ie J dy [U(x, y)A(y, x) - A(x, y)U(y, x)} = 0. (6.9) 

c 

Note that the integral term of this identity is independent of the contour placement of 
the x variable due to discontinuity relation (3.8) and, therefore, it is only a function of 
the physical value of x. 

By means of equations of motion (3.11), (3.19) and (3.20), the divergence of the Noether 
current of eq. (3.27) is seen to vanish 

id,f = e [J*(x)<f>(x) - J(x)<f>*(x)\ +ie J dy [U(x, y)A(y, x) - A(x, y)U(y, x)] = 0, (6.10) 

c 

according to eq. (6.9). Thus, we have arrived at the current conservation for any $- 
derivable approximation, which is invariant under (6.7). 

Similarly, one may derive the relation, resulting from the Lorentz invariance of the $ 
functional, which permits to demonstrate the conservation of the angular momentum. 
However, we do not consider it here, since the angular-momentum conservation is not of 
practical use in kinetics. 

Further invariances generally depend on the properties of the interaction vertices in the 
theory considered. For instance, other invariances result from the property of fixed number 
of field operators attached to the interaction vertices in a theory. This property provides 
relations between densities of the interaction S mt (x) and potential S pot (x) energies. If 
the number 7 of field operators per vertex is fixed, the functional $ is invariant under 
the simultaneous local scaling of the interaction vertices by A(x), and the classical fields 
and propagators by (X(x))^ 1 ^' and (\(x)\(y))~ 1 ^ , respectively. This invariance provides 
relation (2.6) at the level of expectation values. 



7 Thermodynamic Consistency 

In the thermal equilibrium the density matrix is explicitly known, cf. [35], 

^ q exp (-/?£{/*}) 
P Z 
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(7.1) 



where (5 — 1/T is the inverse temperature, and Z is the partition function which is directly 
related to the thermodynamical potential, 



tt = -T\nZ. (7.2) 

Since we deal now with thermodynamics, we have introduced the chemical potential \i in 
the conventional way, i.e. by adding to the Hamiltonian the relevant term 



H{fl}=H-J dVJ(Noether)(z), (7.3) 

where j( Noether ) is the time-component of the charged current eq. (2.10), now with e — 1. 

We can use the same trick as that in the Matsubara technique, i.e. use the fact that the 
equilibrium density matrix formally coincides with evolution operator in the imaginary 
time. In the definition of the W functional (4.2) we explicitly write Tr(F q ... instead of 
(...). Thus, taking into account that r = W at vanishing external sources, we arrive at 
the following form of T functional in equilibrium 



r e V0,0*,A,A,/i} = -iln 



—Trexp 

Zj 



X 



-eq 



T c exp 



-eq 



,(7-4) 



with the integration contour C eq now being the sum of the real-time Schwinger-Keldysh 
contour (see figure 1) and the imaginary-time Matsubara contour, i.e. it starts from an 
initial time to goes to infinity, then back to this initial time and after that, to to — i/?- Taking 
into account the fact that T = for the physical values of 0, 0*, A, and ^-independent 
A, we obtain for the value of the thermodynamic potential (7.2) 



il{(f), (f>* , A, \, /j,} = -ThWTr 





f 










|tt 


exp 


-i J dtH°M 


T c exp 


i J dx\C\ nt {fi} 


)} 






Ceq 









(,7.5) 



where the integral over the real-time section of the contour gives zero. Hence, in eq. (7.5) 
we can make the replacement 



-i/3 



dt... 



dt. 



Ceq 



(7.6) 



Thus, we have arrived at the proper thermodynamic representation of the thermodynamic 
potential originally proposed by Luttinger and Ward [28]. Indeed, since all quantities 
under the integral are analytically continued from the Schwinger-Keldysh contour to the 
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Matsubara contour, Q is determined by the same expression as the T functional (4.7) 
but in terms of the Matsubara Green's functions with the thermodynamic $r functional 
represented by the same set of closed diagrams. Thus, in the momentum representation 
from eq. (7.5) we arrive at 



n 



{0,0*,A,A )/U } = - j dVi^^KT^^ J2 j d 3 x^exp(i^) 

x (- ln[- A(u n , p)} + n(u n , p)A(cu n , pfj + $ T , rj -> 0, (7.7) 



where $t = — iT$, u n = 2irmT, and summation runs over Matsubara frequencies. In 
the standard way (see, e.g., ref. [29]) by converting the w n -sum in eq. (7.7) into the 
energy integral expressed in terms of the real-time quantities of eqs (B.2) and (B.3), this 
thermodynamic potential is also easily expressed in terms of the real-time quantities 



n {0, 0*, A, A, p} = - f d 3 x£°{0, d^} [ d 

J \ Zi J ncut . J 



x———n{e) 



x — 21m In 



ncut. J (27r) Z 

-A R (e + iO, p)] - ReA R T - AReU R ) + $ T , (7.8) 



where 



n(e) = [exp(e/T) - l]" 1 



(7.9) 



is the thermal Bose-Einstein occupation number. Thus, the problem of the thermody- 
namic consistency can immediately be re-addressed from the Schwinger-Keldysh approach 
to the Matsubara one. Within the Matsubara formalism, this problem was considered by 
Baym [27]. He has shown that any ^-derivable approximation to the thermodynamic 
potential is thermodynamically consistent. Hence, we have proved that our $ derivable 
approximations to the T-functional are also thermodynamically consistent. 

The stationary property of the T functional (and, hence, of Q) with respect to variations 
in full Green's functions and classical fields, eq. (4.6), is the key feature of T that provides 
the thermodynamic consistency. It implies that any derivative of the thermodynamic 
potential to any thermodynamic parameter like (3 or \x is then given by accounting only 
the explicit dependence of Q on these parameters, since the implicit dependences through 
A and <fi drop out due to the stationarity property. Therefore ^-derivable approximations 
preserve the corresponding thermodynamic relations as for the exact partition sum, and 
thus provide thermodynamic consistency. 
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8 Conclusion 



Our aim was to develop a regular way of constructing self-consistent approximations to 
quantum transport theory within the formalism of non-equilibrium Green's functions on 
the real-time contour, developed by Schwinger, Kadanoff, Baym and Keldysh [1-3]. 

We have employed functional methods for Green's functions on the real-time contour. 
Our main result is the definition of a generating functional $ which is determined by 
the sum of all closed (i.e. without external points) skeleton diagrams in terms of classical 
fields and full Green's functions on the real-time contour. The main feature of this $ 
functional is that it plays a similar role as the interaction Lagrangian but in the space of 
classical fields and full Green's functions on the contour. This means that all important 
quantities of a system (such as sources of classical fields, self-energies, interaction energy, 
etc.) can be derived by respective variations of the $ functional. Within thermodynamics, 
the $ functional was introduced by Luttinger and Ward [28] and later used by Baym [27]. 
Our treatment extends the definition of the $ functional to any non-equilibrium system. 
Following the thermodynamic treatment of Cornwell et al. [34], we have also extended the 
non-equilibrium <3> functional to the case of non-vanishing classical bosonic fields. This 
last generalization allows us to self-consistently describe the dynamics of both the order 
parameter (the classical field) and fluctuations on equal footing, e.g. in the theory of 
phase-transition phenomena. 

The advantage of the $ functional is that we may formulate various approximations at 
the level of <3>, thus defining so called ^-derivable approximations. In particular, we may 
construct effective theories right at the level of Green's functions and effective vertices. 
These approximations possess some important features: they respect exact conservation 
laws on the level of expectation values (with the Noether values for the conserved quan- 
tities) and have a proper thermodynamic limit. Note that other approximation schemes, 
e.g. at the level of self-energies, far not always possess such properties. 

The question of consistency becomes especially important for a multi-component system, 
where the properties of one species can change due to the presence of interaction with 
the others and vice versa. The "vice versa" is very important and corresponds to the 
principle of actio = re-actio. This implies that the self-energy of one species cannot be 
changed through the interaction with other species without affecting the self-energies of 
the latter ones also. The ^-derivable scheme offers a natural and consistent way to account 
for this principle. Within thermodynamic considerations this has recently been considered 
for the interacting pion - nucleon -delta-resonance system, where the coupling to the delta 
resonance leads to a softening of the pion modes below the resonance mass [36] and for a 
relativistic QED plasma in [37]. 

For the relativistic scheme considered here we argue that a careful construction of con- 
served quantities requires symmetric expressions in terms of A h and A H Green's func- 
tions (A < and A> in the Kadanoff-Baym notation, respectively). This is in contrast to 
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expressions involving only A~ + Green's function, which are often used in the literature. 
These symmetric expressions describe contributions of both particles and anti-particles 
on equal footing, as well as take a proper account of modification of the vacuum polariza- 
tion in the medium. Of course, these symmetric expressions still require a proper vacuum 
renormalization to be done in any actual calculation. 

We have done our consideration on the example of relativistic bosons, since they allow us 
to demonstrate the basic features of this approach — inclusion of both Green's functions 
and classical fields into the scheme — without extra complications. The generalization to 
multi-component systems is straightforward. Here, we expect a consistent description of 
chiral a-, is- condensates together with fluctuations, as an immediate application. The gen- 
eralization to vector mesons and/or fermions, though more technically involved due to the 
resulting tensor structure of the propagators and self energies is also straight forward. The 
formalism can also be generalized to theories with derivative couplings in the interaction 
Lagrangians. While the expressions for the equations of motion and the energy-momentum 
tensor receive additional terms from the derivative couplings, the diagrammatic rules are 
the same as given here. 

The developed scheme of constructing self-consistent approximations provides a suitable 
basis for the derivation of kinetic equations beyond the limitations of the quasi-particle 
approximation. Such generalized transport schemes respect parts of the quantum nature 
of the particles and, in particular, take account of their finite mass-widths. The finite 
mass-width may be either inherent in a particle already from its vacuum properties (e.g., 
resonances) or may be acquired by a stable particle in a dense environment due to fre- 
quent interactions. In the case of nuclear collisions at intermediate (~ 1 GeV/nucleon) 
to ultra-relativistic energies, one encounters mean single-particle energies in the range 
of the typical temperature of T = 50 - 200 MeV. Important resonances, like the delta- 
resonance or the rho- meson, have decay widths beyond ~ 100 MeV, while typical collision 
rates estimated from presently used quasi-particle transport schemes are typically in the 
order of T. These circumstances definitely prevent quasi-particle based transport codes 
from providing reliable results for such collisions. The main steps in the derivation of 
self-consistent and tractable transport equations for particles with finite width will be 
published in a forthcoming paper [38]. 
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APPENDICES 



A Diagram Rules for $, J, and II 

The interaction vertex function V(x) entering the diagram is normalized in the standard 
way, c.f. [32], i.e. with factors n\ relative to C mt (x) for each type of operator occuring with 
multiplicity n in the vertex. E.g., the vertex function simply becomes — iV(xk) = —ig for 
C' mt = -g^/Al (4 identical operators) and for £ int = -g(<p*<j)) 2 / '(2! • 2!) (twice two 
identical operators). The diagrammatic rules to calculate i$, iJ(x) and —ill for a given 
theory are as follows 

(i) Draw all topologically distinct, closed and entirely connected diagrams with N internal 
vertices x±, x 2 , ...xn, where classical field pins and propagator lines saturate the valences 
of all vertices in the diagram, c.f. (5.11) above. Closed diagrams for i$ have no external 
points, while iJ(x) has one external point, and —ill has two external points. For charged 
bosons, pins and propagator lines have an arrows sense, distinguishing <fi from 0^ at 
the vertices, the sense direction pointing towards <fi. 

(ii) For i<3>, iJ{x) and — i£ mt (a;) keep only those diagrams that are two-particle irreducible, 
i.e. which cannot be split into two pieces by cutting two different propagator lines. For 
—ill keep only those diagrams which result from <3> by opening one propagator line. 

(iii) To each line, connecting x\ — > Xk, assign the factor iA(xk,xi). 

(iv) To each pin attached to assign the factor <f>(xk) or <p*(xk) depending on the sense. 

(v) To each vertex Xk assign the vertex factor —iV(xk) as determined by C int (x). 

(vi) Integrate all internal xi,x 2 , ...xn over the contour. 

(vii) Multiply the result by the symmetry factor S, which is calculated as follows 
1/iVc?! factor for every Nq equivalent internal lines, 

l/N^. factor for every classical fields entering each vertex, 
1/2 factor for every self-closed line loop (tad-poles) for real fields. 

(viii) Sum all diagrams. For i$ (contrary to iJ(x) and —ill), put the extra factor \jn\ for 
each diagram, where ri\ counts the number of vertices in the diagram. 

In many cases like in transport treatments, it is advantages to consider the diagrams 
decomposed into the two contour sections at each vertex, e.g., to calculate quantities like 
II h and Il + ~ self-energies in terms of exact Green's functions. Therefore, "physical" -time 
diagrammatic rules in the matrix scheme are also required. Here we present only those 
rules which differ from the above ones on the real-time contour, bearing in mind that all 
other rules remain valid: 

(iii') To each internal vertex Xk first assign a sign %\. e {+, — } defining the contour placement 
x l £. To each line, connecting x\ l — > x l £, assign the factor iA tkt '(xk, x{), ik,h £ {+, — }• 

(vi') For all internal points integrate all xi, x 2 , ...xn over the real-time axis and space, for 
each internal " + " vertex multiply by (—1) and finally sum over all internal contour 
placements ii, i 2 , ■■■In {ik € {+, — }). 
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B Equilibrium Relations 



For completeness of the thermodynamic consideration, we explicitly present here equilib- 
rium relations between quantities on the real-time contour. Basically, they follow from 
the Kubo-Martin-Schwinger condition [39] 

A- + (p) = A+-(p)e-c/ T , U- + (p) = U + -(p)e-^ T , (B.l) 

where e = p u U u — fi with U" and p, being a global 4-velocity of the system and a chemical 
potential related to the charge, respectively. All the Green's functions can be expressed 
through retarded or advanced Green's functions: 



/ [1 + n{e)\ A R (p) - n(e)A A (p) -m{e)A{p) \ 

(A^ip)) = 

\ -i [1 + n(e)} A(p) -[l + n(e)]A A (p)+n(e)A R (p)J 

i,j e {+, — }, and the self-energies take a similar form 

( U. R {p) -\n{e)Y{p) -ira(e)r(p) N 
v -i [1 + n{e)\ T{p) -U A (p) - m(e)r(p) ) 



(,B.2) 



(U^(p)) = 



(B.3) 



Here n(e) is the thermal Bose-Einstein occupation number defined in eq. (7.9). 
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